Terahertz Quantum Hall Effect in a Topological Insulator 
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Using THz spectroscopy in external magnetic fields we investigate the low-temperature charge 
dynamics of strained HgTe, a three dimensional topological insulator. From the Faraday rotation 
angle and ellipticity a complete characterization of the charge carriers is obtained, including the 2D 
density, the scattering rate and the Fermi velocity. The obtained value of the Fermi velocity provides 
further evidence for the Dirac character of the carriers in the sample. In resonator experiments, we 
observe quantum Hall oscillations at THz frequencies. The 2D density estimated from the period of 
these oscillations agrees well with direct transport experiments on the topological surface state. Our 
findings open new avenues for the studies of the finite-frequency quantum Hall effect in topological 
insulators. 



Three dimensional topological insulators 0, Q have at- 
tracted much interest recently, as they exhibit a number 
of unusual and non-trivial properties, such as protected 
conducting states on the surfaces of the sample. Un- 
usual electrodynamics, such as a universal Faraday effect 
and an anomalous Kerr rotation have been predicted 0- 
6j for these surface states, their observation is still out- 
standing. We showed recently that strained HgTe, where 
the strain lifts the light-hole-heavy-hole degeneracy that 
normally is present in bulk HgTe, is a very promising 3D 
topological insulator [?J • This is because at low tempera- 
tures parasitic effects due to bulk carriers are practically 
absent. In static transport experiments a strained 70 
nm thick HgTe layer Q exhibits a quantum Hall effect 
(QHE), yielding direct evidence that the charge carriers 
in these layers are confined to the topological two dimen- 
sional (2D) surface states of the material. These findings 
are further corroborated by recent Faraday rotation data 
Q in a similar layer, which have been obtained using a 
terahertz time-domain technique. 

In this work, we present the results of low tempera- 
ture terahertz Faraday cw transmission experiments on 
another strained HgTe film. The carrier density, Fermi 
velocity and the scattering rate can be reliably deter- 
mined from these data. In particular, we obtain the 
Fermi velocity vp = 0.52 • 10 6 m/s, which is in excellent 
agreement with the Faraday rotation experiments 0] and 
the dc Shubnikov-de Haas measurements Q on 70-nm- 
thick strained HgTe films as well as with band-structure 
calculations for the surface states in 3D topological insu- 
lators (see e.g. Ref. Q). In the same sample we observe 
quantum Hall-induced oscillations at terahertz frequen- 
cies, providing further evidence for the 2D character of 
the conductivity. In the case of topological insulators, no 
finite frequency QHE has been reported up to now. 

The sample studied in this work is a coherently 
strained 52-nm-thick nominally undoped HgTe layer, 



grown by molecular beam epitaxy on an insulating CdTe 
substrate (lo| . Transmittance experiments at terahertz 
frequencies (100 GHz < v < 800 GHz) have been 
carried out in a Mach-Zehnder interferometer arrangc- 
[l2j which allows measurement of the ampli- 
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tude and phase shift of the electromagnetic radiation in 
a geometry with controlled polarization. Using wire grid 
polarizers, the complex transmission coefficient can be 
measured both in parallel and crossed polarizers geom- 
etry Static magnetic fields, up to 8 Tesla, have been 
applied to the sample using a split-coil superconducting 
magnet. 

To interpret the experimental data we use the ac con- 
ductivity tensor cr(w) obtained in the classical (Drudc) 
limit from the Kubo conductivity of topological surface 
states (see e.g. Ref. Q). The diagonal, cj xx {uj), and 
Hall, a xy (uj), components of the conductivity tensor as 
functions of THz frequency ui can be written as: 



<7xx(u) = (Jyy(uj) = 

Vxy{u) = -<T yx (uj) ■ 



1 — IUJT 



(1 - {.cur) 2 + (O c r) 2 



(1 - iujT) 2 + (n cT ) 2 



<?v , (1) 
cr . (2) 



Here, fl c = eBvp/hkp is the cyclotron frequency, oq is 
the dc conductivity, B is the magnetic field, Vp, kp, e, 
and r arc the Fermi velocity, Fermi wave-number, charge, 
and scattering time of the carriers, respectively For the 
Dirac spin-helical surface states the Fermi wave-number 
depends on the 2D carrier density, n^B 1 through relation 
kp = \fAnri2D, with no spin degeneracy. 

The transmission spectra can then be calculated using 
a transfer matrix formalism [uj - froT ] which takes multiple 
reflection within the substrate into account. The electro- 
dynamic properties of the CdTe substrate have been ob- 
tained in a separate experiment on a bare substrate. Fur- 
ther details of the fitting procedure can be found in the 
Supplementary information to Ref. [15j. Neglecting any 
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FIG. 1: Magnetic field dependence of the transmission in 
strained HgTe. (a-c) Transmission amplitude in parallel po- 
larizers (t p ) geometry, showing cyclotron resonance at the 
positions indicated by the arrows. The frequency of the ex- 
periments is indicated in the panels. The inset shows the 
frequency dependent transmittance in zero external magnetic 
field, \t p (B — 0)j 2 . Symbols: experiment, solid lines: simul- 
taneous fit of all data with the Drude model as described in 
the text. 



substrate effects, the complex transmission coefficients in 
parallel (t p ) and crossed (t c ) polarizers geometry can be 
written as: 



in 
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Here T, xx and T, xy are effective dimensionless 2D con- 
ductivities, defined dZ 
with the HgTe film thickness d = 52 nm and the vac- 
uum impedance Zq « 377 O. In order to self-consistently 
obtain the parameters of the quasiparticles, the field- 
dependent complex transmission t p (B) and t c (B) for 
v =0.17 THz, 0.35 THz and 0.75 THz and the zero-field 
transmittance spectra |£ c (w)| 2 have been fitted simulta- 
neously. 

The inset in Fig. [1] shows the transmittance spectrum 
of the HgTe film at zero magnetic field. The character- 
istic oscillations in the spectrum, with a period of about 
58 GHz, are due to Fabry-Perot type interferences within 
the CdTe substrate. The absolute transmittance in the 
interference maxima is close to 95%, which reflects the 
low effective conductance of our HgTe film, T, xx <C 1. At 
low frequencies, the maximum transmittance decreases 
and approaches |i p | 2 ~ 0.7 in the zero frequency limit. 
Such a behavior is typical for Drude carriers with a scat- 
tering rate in the frequency region of the experiment. 
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FIG. 2: Complex Faraday angle 6 + in in HgTe. Bottom 
panels: Faraday rotation, top panels: ellipticity for the same 
frequencies as in Fig. [TJ The inset sketches the definitions 
of the Faraday rotation 9 and ellipticity 77. Symbols: experi- 
ment, solid lines: simultaneous fit of all data with the Drude 
model as described in the text. Angular units are radians. 



Indeed, the solid line in the transmission spectra repre- 
sents a Drude fit with the parameters given in the first 
row of Tab. H 

From the fits we obtain the Fermi velocity vp = 
0.52 • 10 6 m/s. This value is very close both to vf = 
(0.51-^0.58) • 10 6 m/s as determined in the Faraday rota- 
tion experiments on a 70-nm-thick strained HgTe film 
[H and to vf — 0.42 • 10 6 m/s as extracted from dc 
Shubnikov-de Haas measurements on a patterned 70-nm- 
thick strained HgTe layer Q. The obtained value of the 
Fermi velocity is also in very good agreement with the 
band-structure-theory result vp = 0.51 • 10 6 m/s for the 
linear (Dirac) part of the surface-state spectrum in topo- 
logical insulators (see e.g. Ref. Q). As an additional 
check of the 2D surface carrier dynamics in our sam- 
ple, we have analyzed the terahertz transmission data of 
Ref. 15( for a 70-nm-thick strained HgTe film at high 
temperature T = 200 K and for a bulk (1000-nm-thick) 
unstrained HgTe sample. In both cases, the electrody- 
namics is governed by massive bulk carriers, for which 
the values of Vp turn out to be much larger than the 
Dirac surface-state velocity, i.e., vf ~ 0.5 • 10 6 m/s (Tab. 



Figure [T] shows the magnetic field dependent trans- 
mittance of the HgTe film in Faraday geometry and for 
parallel orientation of polarizer and analyzer. According 
to Eq. ([3]), the transmittance in parallel polarizers (t p ) 
depends mainly on Y> xx . For all three frequencies two 
clear minima in the transmitted signal are observed in 
the range below ±1 T. The minima in \t p \ roughly corre- 
spond to the cyclotron resonance energy and scale with 
magnetic field. This may be understood taking into ac- 
count that in our case E < 1 and Eqs. (|3I4[) simplify 
to: 
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TABLE I: Drude parameters of the charge carriers in HgTe in strained and unstrained films. The data on 70 nm and 1000 nm 
film were partly given in Ref. [l5l |. 



Sample 


r(K) 


n 2 D(cm 2 ) 


vf(ths ) 


1/27TT (GHz) 


G 2D =ao-d (IT 1 ) 


52 nm (strained) [this work] 


2 


1.08 ■ 10 11 


0.52 ■ 10 b 


250 


7.6- 10~ 4 


70 nm (strained) [15] 


4 


4.8 • 10 iu 


0.38 • 10 b 


210 


4.3 ■ 10" 4 




200 


1.5 ■ 10 12 


1.63 • 10 6 


360 


5.3 ■ 10" 3 


1000 nm (unstrained) [15] 


3 


4.2 • 10 11 


0.99 ■ 10 b 


240 


2.8 • 10" 3 




200 


4.9 ■ 10 13 


9.36 • 10 6 


360 


1.9 ■ lO^ 1 



In the limit ujt 3> 1, Eq. ([T]) may be approximated by 



(0) 



which leads to a resonance like feature for O c = uj. Thus, 
the positions and widths of the minima in Fig. [T] are 
directly connected with the parameter vfI^f and the 
scattering rate r _1 of the charge carriers. 

Figure [2] shows the complex Faraday angle 9 + irj as 
obtained at the same frequencies as in Fig. [TJ The po- 
larization rotation 9 and the ellipticity r\ are obtained 
from the transmission data using: 



tan(20) = 25R( X )/(lHx| 2 ) 
sin(2 ?? ) = 23( X )/(1 + | X | 2 ) 



(7) 

(8) 



Here x = t c /t p and the definitions of 9 + irj are shown 
graphically in the inset to Fig. [5J A direct interpretation 
of the complex Faraday angle is in general not possible 
because of the interplay of a xx and a xy in the data. 

In the low frequency limit, lot <C 1 Eq. ([2]) simplifies 
to the static result a xy = ft c Ta n /(l + (f2 c T) 2 ). The last 
expression has a maximum at fl c (B) = r , which leads 
to maxima in t c and 9 at about the same field value. 
Therefore, the Faraday angle provides a direct and an 
independent way of obtaining the scattering rate 1/r. 
The solid line in Fig. [5] are the fits which have been 
done simultaneously for all results presented above. In 
total, the parameters of the charge carriers have been 
obtained by simultaneously fitting ten data sets. The 
quite reasonable fit of all results proves that a single type 
of charge carriers dominates the electrodynamics in the 
range of frequencies and magnetic fields used in these 
experiments. 

Very solid evidence for the two dimensional charac- 
ter of the carriers probed in the Faraday rotation ex- 
periments would be the observation of quantum Hall 
plateaus, similar to the observation of the QHE in Q. 
However, the accuracy of the experiments shown above 
does not allow to observe the QHE. In order to solve 
this problem, we have performed further Faraday trans- 
mission experiments on the same sample, now using a 
resonator geometry as shown in the inset of Fig. [3] 

In these experiments, the sample is placed in the 
middle of a Fabry-Perot resonator defined by metallic 
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FIG. 3: Faraday rotation in HgTe within resonator geome- 
try, (a) - Ellipticity, (b) - Faraday angle. Symbols - ex- 
periment, lines - fits according to Eqs. (1-4). Upper inset 
shows the experimental geometry within a Copper meshes 
resonator. Lower inset shows a magnified view of the Fara- 
day angle demonstrating QHE oscillations. 



meshes. We have utilized Cu meshes with a 200 /iin 
period. The distance between adjacent maxima of the 
resonator is ~ 51 GHz. In the frequency range between 
100 and 200 GHz the quality factor of the loaded res- 
onator is about Q ~ 10. This indicates that, effectively, 
the radiation passes about ten times through the sample 
before reaching the detector, which effectively increases 
the sensitivity to fine details by roughly the same value. 
As shown in Fig. [3J in the resonator experiments the 
field dependence of the Faraday rotation and the ellip- 
ticity appears qualitatively similar to that in Fig. [2] 
An exact calculation of the complex transmission coef- 
ficients within a resonator is complicated because of the 
increased number of parameters. Therefore, in this case 
we utilize the simple equations Eqs. (JT])-(|4j) which neglect 
the effect of the substrate and the resonator completely. 
Nevertheless, as clearly seen in Fig. 02 the fits based on 
the simplified expressions reproduce the experimental re- 
sults reasonably well. Fitting of the signals for parallel 
and crossed polarizers yields within experimental accu- 
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FIG. 4: Terahertz quantum Hall effect in HgTe. (a) 
Two dimensional conductance: G xx , (b) derivative of G^y 
(dGxy/dB- 1 ). The data have been obtained within a res- 
onator geometry and are plotted as a function of inverse mag- 
netic field. The experimental data are shown as solid lines for 
frequencies as indicated. Dashes in the bottom panel marks 
the minima for negative magnetic fields, (c) Numbered posi- 
tions of the minima in G xx and in the derivative of G xy for 
0.14 THz and 0.19 THz. Straight lines yield interpolation to 
the origin. 



racy the same parameters as in the experiments without 
a resonator. The only parameter which differs from the 
results without a resonator is the absolute value of the 
conductivity. This is of course expected, and results from 
multiple transmission in the resonator and the influence 
of the substrate. 

The main advantage of the resonator experiments is 
a higher sensitivity to the details of the field-dependent 
transmission. In addition to an overall field dependence 
similar to that in Figs. [1] and [2] a tiny modulation of 
the signal can now be observed. To convert this modu- 
lation to a conventional presentation, we have inverted 
the transmittance curves into the 2D conductivity, using 
Eqs. ([3]and|4]). Because the absolute transmittance is 
not well-defined in the resonator experiments, we have 
scaled the absolute 2D conductance to agree with the 
data without a resonator. The final results expressed in 
form of the effective 2D conductance G XXiXy = ^ xx , xy /Zo 
are shown in Fig. |4l 

Fig. 0k shows the real part of the two dimensional 
conductance G xx as a function of inverse magnetic field. 
Clear oscillations in the conductance can be observed 
in this presentation. In general, the phenomenology of 
the QHE at terahertz frequencies is not well understood 



la . Il7| . Existing experiments are generally limited to 



frequencies below 100 GHz and they are analyzed using 
scaling exponents flil [l8j]. In the resonator experiments, 
the field dependent oscillations can be observed both with 
parallel and crossed polarizers. Contrary to G xx , the 
off-diagonal conductance G xy shows a substantial field 
dependence even in high magnetic fields. Therefore, no 



clear QHE signal can be directly detected in G xy . In or- 
der to extract the QHE information from these data, we 
have plotted the derivative of the G xy as a function of 
an inverse magnetic field (dGxy/dB -1 ) in Fig. HJd. The 
derivative has the advantage of being insensitive to any 
residual slowly varying signals, and, importantly, the ex- 
pected steps in G xy are transformed into the minima of 
the derivative. Finally, in order to analyze the quantum 
Hall effect, both the minima in G xx and in the derivative 
of G xy have been taken into account. In Fig. Hk,b the 
results at finite frequencies are compared with dc QHE 
on the same sample. The periodicity of the oscillations in 
the dc experiments is slightly different because of differ- 
ent carrier concentration at the sample surface, induced 
by exposure to photoresist and the presence of ohmic 
contacts. 

The main results of the QHE experiments are repre- 
sented in Fig. HJ: demonstrating an approximate equidis- 
tant positioning of all minima (labeled by number N) in 
inverse magnetic fields B~ x with the period of AB^ 1 = 
0.18 T _1 . This periodicity reflects the dependence of the 
number of the occupied Landau levels on B . In a to- 
tal, we detect the oscillations up to index number ±10; 
also the first oscillations with the Landau level index ±1 
are clearly observed in the data. From the periodicity of 
these oscillations the effective 2D carrier density can be 
estimated according to the free electron expression n 2 D = 
e/(ft,A_B _1 ) ~ 1.4 • 10 11 cm -2 . This value agrees reason- 
ably well with the density n2D = 1.08-10 11 cm -2 obtained 
directly from fitting the transmittance and the Faraday 
rotation on the basis of the Drude model (Tab.|T|. There- 
fore, we may conclude that charge carriers which are re- 
sponsible for the terahertz electrodynamics at low tem- 
peratures reveal 2D behavior. To further characterize the 
electron system in our sample we extrapolated the de- 
pendence N(B~ 1 ) to the origin (see straight lines in Fig. 
2J:), which corresponds to the limit of very strong mag- 
netic fields. At the origin we find a finite value JV k 1/2 
instead of N = as would be the case for the conven- 
tional QHE. Previously, similar extrapolated values were 
reported for graphene (see e.g. Ref. [l9|), zero-gap HgTe 
quantum wells [20( and strained 70 nm-thick HgTe films 
[Jj, i.e. for materilas with 2D Dirac-like charge carriers 
encoding a nonzero Berry phase. We therefore believe 
that our teraherz QHE also indicates the 2D Dirac-like 
behavior. 

In conclusion, we have analyzed the terahertz Fara- 
day rotation in a strained HgTe film. From these data 
all relevant parameters of the charge carriers can be ob- 
tained. In addition, terahertz quantum Hall effect oscil- 
lations have been observed within the same experiment, 
which proved the two-dimensional character of the con- 
ductivity. 
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